Marginal densities of the "true" self-repelling motion 



Laure Dumaz and Balint Toth 



February 21, 2012 



Abstract 

Let X{t) be the true self-repelling motion (TSRM) constructed in |20j . L{t,x) 
its occupation time density (local time) and H{t) := L{t,X{t)) the height of the 
local time profile at the actual position of the motion. The joint distribution of 
H{t)) was identified in |18) in somewhat implicit terms. Now we give explicit 
formulas for the densities of the marginal distributions of X{t) and H{t). The 
distribution of X{t) has a particularly surprising shape (see Picture ([1]) below): It 
has a sharp local minimum with discontinuous derivative at 0. As a consequence 
we also obtain a precise version of the large deviation estimate of [1]. 
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Figure 1: Density of X[l) 



Figure 2: Density oi H(l) 



1 Introduction and main results 



1.1 Introduction 

In the present paper, we study some marginal densities of a self-interacting one dimen- 
sional process called the true self- repelling motion (TSRM), defined in [20]. The true 
self-repelling motion is a continuous real- valued process {X{t), t > 0) that is locally 
self-interacting with its past occupation-time: It is pushed away from the places it has 
visited the most in its past. This process has a completely singular behavior compared 
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to diffusions: For example, it has a finite variation of order 3/2, which contrasts with the 
finite quadratic variation of semi-martingales. 

A crucial property of the TSRM is the existence of a continuous density for its 
occupation-time measure, denoted L{t,x) and named again "local time" by analogy with 
semi-martingales (see Theorem 4.2 of pO]). This enables to prove that the TSRM is 
driven by the negative gradient of its own occupation time density L{t,x). Written 
formally: 

dtX{t) = -d,L{t, dtL{t, x) = 6{x - X{t)) (1) 

X(0) = 0, L{0,x) = 0. (2) 

The driving mechanism ([1]) is to be considered in a properly regularized form. Making 
proper mathematical sense of this mechanism is the main content of |20) . 
The local time at the current position of the TSRM is denoted by 

H{t) ■.= L{t,X{t)) 

and is called the height process. The processes X{t), the local time L{t, x) and the height 
H{t) scale as follows: 

X{at) ~ a2/3x(t), L{at, a^/^x) ~ a^^^L{t, x), H{at) ~ a^/^H{t). 

We will examine in this paper the distributions of the marginals X{t) and H[t) at the 
fixed time t. It turns out that it is more convenient to work at first with an exponential 
random time of parameter s independent of the process instead of with a fixed time. 
With Feynman-Kac formulas applied for the Brownian motion, it is possible to deduce 
explicit expressions for its marginals (in terms of Airy function, confluent hypergeometric 
functions and Mittag-Leffler density). It leads to the Pictures [T] and |2] for t = 1. We 
found the shape of the of the density of the position X{t) rather unusual and surprizing. 
Indeed, it has a sharp wedge-like local minimum at a; = with discontinuous derivative, 
and global maxima away from zero. It means that at any positive time, the process still 
remembers its starting point and it is strongly pushed away from it. These results are 
also confirmed by numerical simulations (see Figures [31 and HI). 

1.2 Review of background and notations 

The true self- avoiding walk (TSAW) is a nearest neighbor random walk with long memory 
on 1/ pushed by the negative gradient of its own occupation time measure (local time). 
For historical background see j2], jT3], |T3] or the introductions to [18j or [20j. In the 
present paper we are interested only in the one-dimensional case. 

A discrete time version considered and analyzed in [18] is the following. Let n ^ 
S{n) G Zi be a nearest neighbor walk on Z. The occupation time on lattice bonds, lattice 
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sites, and the negative (discrete) gradient of bond occupation time are in turn: 

i{n,j + 1/2) := #{0 < m < n : {S{m) + S{m + l))/2 = j + 1/2}, 

i{n,j) := 1/2) +£(n,j + 1/2)) /2, 

6{n,j) := i{n,j -1/2) ~i{n, J + 1/2), 

where j G Z label lattice sites, J ± | G ^ + i label lattice bonds and n G is the 
(discrete) time of the process. 

The TSAW (with bond repulsion) is defined by the law: 

P {S{n + 1) = J ± 1 I S{n) = j) = "".[tT'Pi 

I w{6{n,j)) +w{-6{n,j)) 

where if : Z — )■ (0, oo) is a fixed rate function which is assumed nondecreasing and non- 
constant. In jTS] limit theorem was proved for the distribution of the displacement of 
the TSAW, with time-to-the-two-thirds scaling. Later, the convergence was established 
for the whole process [20| [12]. In order to expose this and also prepare our present work 
we need some notation and preliminary concepts. 

Let M 9 y I— )■ B{y) G M+ be a two sided Brownian motion, starting from level h > 0, and 
let u' be the first hitting time of zero by the backward Brownian motion, the first one 
after time x > of the forward Brownian motion: 

u' := sup{y < : B{y) = 0}, (3) 

CO, := mf{y > x : B{y) = 0}, x> 0, (4) 

and denote: 

T,= / \B{t)\dt. (5) 



The random variable is a.s. positive and finite and has an absolutely continuous 
distribution with density 

g{t;x,h) := dtP{T^ < t \ B{0) = h), te R+, x eR+, he M+. 

The Laplace transform of this density with respect to the variable t is 

POO 

g{s;x,h):=s e''^g{t;x,h)dt = sE{e~'^^ \ B{0) = h). (6) 
Jo 

We extend g and ^ to x G M as even functions: g{t;x,h) = g{t; \x\,h) and g{s;x,h) = 
g{s] \x\,h) for all x G M, > 0, s > and t > 0. Let us introduce the following marginals: 

/•oo fCO 

gi{t;x):= / g{t;x,h)dh, g2{t]h) := I g{t;x,h)dx, 



gi{s;x):= / g{s;x,h)dh, g2{s;h):= / g{s;x,h)dx. 



oo 
oo 
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Brownian scaling (-B(a-) | B{0) = y/ah^ ~ (^y/aB(-) | -B(O) = /i) implies the following 
scaling relations: 

g{t; X, h) = t-V(r2/3x, t-^/^'h), Qi{t- x) = t-^/^Ui{t-^/^x), Q2{t; h) = t'^'^^ir^''^, 

(7) 

q{s] X, h) = su{s^/^x, s^/'^h), gi{s; x) = s^/^i)i{s^/^x), g2{s; h) = s^l^U2{s^'^h), 
where 

v{x,h) := Q{l]X,h), vi{x) := Qi{l]x), U2{h) := g2{l; h) , 

i'{x,h) := g{l;x,h), z>i(a;) := x), U2{h) := Q2i^;h). 



We recall the main results of [IB] and [20] and more recent large deviation estimates 
from [1| which are necessary background material for the present paper: 

• Theorem 2 of [18j: {x,h) i— )■ v{x,h) and {x,h) ^ j){x,h) are probability densities 
on M X M+: 

/oo /»oo fOO /»oo 

/ z/(a;, h)dhdx = 1=1 / i>(x, h)dhdx. (8) 
oo Jo J -oo Jo 

• Theorem 3 of |18] : Let s > be fixed and 6n a sequence of geometrically distributed 
stopping times with E = n/s, independent of the TSAW S{-). The following 
limit theorem holds: 

where the density of the joint distribution of (X^, Hg^ is g{s] x, h). The constant a 
in the norming factors on the left hand side depend only on the rate function w{-) 
and is given by formula (1.23) in [18J. 

Actually, in [TH] only the case 'w{z) = e^^^ was considered and only the limit 
theorem for the displacement was stated explicitly. However, the proofs implicitly 
contain these extensions. 

• Remark after Theorem 3 of [18j: Assume that the sequence of random variables 
{n-^/^S{n),n-^/^(.{n, S{n))) G M x M+ is tight. Let t > be fixed. Then 

{{an)-'/'S{[nt]),{an)-^/H{[ntlS{[nt\))) {X{t),H{t)), 

where the density of the joint distribution of {X(t), H(t)) is g(t; x, h). 

See also [19] for similar results for a continuous time version with on-site-repulsion. 
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• Finally, in |[2U] the true self-repelling motion (TSRM) introduced in the first para- 
graph of this paper, 

t ^ X{t) G M 

was constructed. This is the scaling limit of the properly rescaled TSAW, t i— )■ 
(an)'"^^^ S {[nt]) , as n — )■ oo. It implies that the density of the joint distribution of 
the position and height processes {X{t),H{t)) is exactly 

dlf^P {X{t) < X, H{t) <h) = g{t; x, h). 

Nevertheless, it is interesting to note that it can be proved directly in the continuous 
setting thanks to the local time properties established in [20]. Indeed, TSRM's 
construction immediately gives a strong "Ray-Knight" Theorem concerning its local 
times at certain random times (see Theorem 4.3 in [20]). Here is a weaker version: 
For every [x, /i) G M x ]R+, let us denote by r^^h the first hitting time of [x, h) by the 
process {X[t), H{t)). Then the local time at time r^^h has the distribution of a two 
sided Brownian motion starting at time x at the level h reflected above the x-line 
between and x and absorbed by the x-line outside this interval. In particular, it 
implies that T^^h has the distribution of conditionally on -B(O) = h. Therefore, 
with the measure preserving change of variable t = T^^h, we can write for every 
bounded Borel function / on R x R+, 

/+00 ;'+oo 
/ fix,h)E{e-'^^ \B{0) = h)dhdx. 
■oo Jo 

The scaling property readily gives that the density of {X{t), H{t)) at (x, h) G ]RxR+ 
is equal to the density of T^^h at time t i.e. g(t; x, h). 

• By choosing the initial conditions 

X(0)=0, L(0,x) = S(x), 

where x H- B{x) is a two sided Id Brownian motion, rather than ([2]), we obtain 
a version of the TSRM with stationary increments. Denote this process Xst{t), its 
local time Lst{t,x), and the height process Hst{t) := Lst{t, Xst{t)). The scaling of 
the stationary version is similar to that of the TSRM starting from initial conditions 

Xst(at) ~ a'/=^Xst(t), Lst(at, a^/^x) ~ a^^^L.tit, x), H.tiat) ~ a^/^H,tit), 
From the construction of TSRM it follows directly that 

X,t(l)~2-V3x(l). (9) 
But no distributional similarities of this sort hold for the Id marginal ifst(l)- 
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• In |3], using TSRM's construction, asymptotics for the tails of X{1) and -f^(l) are 
found. It is shown that: 

P (X(l) >x)= exp (-(45;V)/27 + 0(ln(x))) , (10) 
P {H{1) >h) = exp {-i8h^)/9 + 0(ln(/i))) , 

where 6[ is the absolute value of the first negative zero of the derivative of a renor- 
malized Airy function. It is precisely defined in paragraph 13.21 

In the stationary case, the tails of the distributions are also computed. In this 
case, due to (Q, the constant in front of the is multiplied by 2 for Xst(l). The 
height -ffst(l) behaves differently as the initial local time helps to have large values 
of |iyst(l)|- It is shown that there exists rj > such that for every large h: 

exp (-l/r/Zi^/^) < P(±iy,t(l) > h)< exp {-?] 0^) . (11) 



1.3 Main results 

In the present paper, our main goal is to identify the marginal distributions, z/i(x) and 
i'2{h) more explicitly than it was done in [T8]. Our main results are collected in the 
following theorem: 

Theorem 1. (i) The second marginal densities C'2{h) and V2{h) are 

P2{h) = -2u{h)u\h) = -dh [u") (/i), (12) 

where u : IR+ — )■ M the (rescaled and normed) Airy function, as defined in subsection 
\3.S\ (142|) . ( |43|) ■ and f/(l/2,4/3; ■) is the confluent hypergeometric function of the second 
kind (or Tricomi function) , given by the integral representation (132]) in subsection \3.4 - 

(a) The first marginal densities Pi (a;) and z/i(x) are 

oo ^, 

0,{x) = Y,Pk^exp{-5i\x\), (14) 
fc=i 

oo 

i^iix) = J2p,^f2d6M), (15) 

k=l 

where f2/3 is the Mittag-Leffler density function as defined in subsection the scaling 
factors are the zeros of the derivative of the (rescaled and normed) Airy function u{-), 
and the coefficients pk of the convex combinations are those given in subsection \3.S\ (HHj) . 



6 



It is worth to compare the computed graphs of these density functions with histograms 
of the empirical distributions on X{1) and H{1) (see Figures [31 and|H). 




Figure 3: The histogram represents a numerical simulation of the TSRM at time one 
X(l) with a sample of 10^ independent realizations. Each realization is computed using 
the discrete "toy-model" introduced in pO| with a step of 2 ■ 10~^. The plain blue line is 
the graph of the density of X{1). 




Figure 4: The histogram represents a numerical simulation of the height process at time 
one H{1) (sample of 10^ independent realizations using discrete "toy-model" with a step 
of 2 ■ 10~^). The plain blue line is the graph of the density of H{1). 

The moments of the marginal distributions collected in Theorem [1] are as follows 
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As a direct consequence of f ll3p and f llSp we also obtain the precise tail asymptotics 
of i'2{h) and as h ^ 1, respectively, as |x| ^ 1 

Corollary 1. 

lim h-HogP{H{l) > h) 

h—^oo 

lim a;"^logP {X{1) > x) 

x—>oo 

In view of the tail asymptotics fllSp also implies 

lim log P (Xst(l) > x) 

Note that this is matching with the large deviation results of [4j. More precision in the 
asymptotics can be derived from the expression of the densities. 

The structure of further parts of this note is as follows: In section |2] we present 
the proofs of the main statements. Section [3] is an Appendix which contains classical 
ingredients: the Feynman-Kac formulas used, and collections of facts about Airy function, 
Mittag-LefHer distributions and confluent hypergeometric functions. We do not prove the 
statements collected in the Appendix but give precise reference for all. 

2 Proofs 

2.1 Preliminaries 

Recall the definitions of T^, and u' {x > 0) in II. 2^ (jlj), and We write 
where 

S:= \Biy)\dy, := \B{y)\dy, := / \B{y)\dy, 

Jul' Jo Jx 

(see Picture [51) 

The main ingredients of our investigations are the following three functions: 

u{h) ■.= E{e^^ \ B{0) = h), h>0, (19) 

^{x, h) := E{e-^--^- \ B{0) = h), x>0, h>0, (20) 

w{x) := if{x, 0) = E(e"^-"^' | B{0) = O) x > 0. 
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(0,0) 



Figure 5: Representation of = S + + 



We also define the Laplace transforms (in the variable x > 0). 



(p{X,h) := / e ^''(p{x,h)dx, 



(21) 



w{\) := [ 



oo 



e w{x)dx. 



Jo 



Remark. Note that this is a different type of Laplace transform than the one we used 
for the time variable. That is why we use different notations (hat for the time transform, 
tilde for this last one). 

In the following Proposition we collect the basic starting identities. All these fol- 
low from straightforward applications of infinitesimal conditioning and/or Feynman-Kac 
formula. 

Proposition 1. (i) The function u : R+ — )■ [0, 1] defined in (fT9|) is the unique hounded 
solution of the boundary value problem 



That is: It is exactly the normalized Airy function defined and treated in subsection \3.^ 
restricted to h & [0, oo). 

(a) The function if : ]R_|_ x ]R+ — )• [0, 1] defined in (120]) is the unique bounded solution of 
the parabolic Cauchy problem 



)lu{h) = 2hu{h) 



m(0) = L 



(22) 



d^(p{x, h) = -dl(f{x, h) - hip{x, h), 



(23) 



with initial and boundary conditions 



IC: <^{0,h) = u{h) 



BC : dh<^{x, 0) = 0, for x > 0. 



(24) 



9 



(in) The function : R_|_ x — )• [0, 1] is the unique bounded solution of the parabolic 
Cauchy problem 

d^v{x, h) = [u{hfd^ [uihyvix, h))) , (25) 

with initial and boundary conditions 

IC: u{0,h)=u{hf, BC: u{x,0)=w{x) for x > 0. (26) 
Proof of Proposition Ul 

(i) (122|) follows from usual "infinitesimal conditioning". 

(ii) By strong Markov property of the Brownian motion 

ip{x,h) = E(e-/o'l^(^)l'^%(|5(x)|) I B{0) = h), 

and thus, we apply the Feynman-Kac formula from subsection 13. 1[ with V{h) = \h\ and 
f{h) = u{\h\). Since for all x > 0, /z. i— )■ ip{x,h) is even and smooth, the boundary 
condition in follows. 

(iii) Note that 

z>(x, h) = u{h)(p{x, h), (27) 

and use ((22]) , ([23]) and 1^. 

□ 



2.2 Computation of z>2 

Since 



oo 



z>2(/i) = 2 / j){x,h)dx, 



from f l25]) and f l26]) we readily get 

(M(/i)2(M(/i)-2z>2(/i))y = -4m(/i)2. (28) 

Using the following identities (that directly come from integration by parts and the 
differential equation of u): 

h 

2 J,, u„.ru\2 „,//L\2/o I „.Un\'2 







uixYdx = hu{hy - u\hY/2 + u'(0)72, 

u{x)-'u'{xfdx = -u'{h)/u{h) + w'(0) + h\ 
we deduce that the general solution of (128]) is 

-2u{h)u'{h) + Cim(/i)2 + C2u{hf / u(x)"^rfx- (29) 
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From the asymptotics (1451) it follows that 



/ u{x) "^dx h as /i — )■ oo. 

Jo 

Since, by ([8]), 

POD 

/ U2ih)dh = 1, (30) 
Jo 

we must choose C2 = 0. (Otherwise 1/2 wouldn't be integrable.) Finally, due to (1301) 
again, Ci = 0, too. The proof of (II 2p is completed. 

Remark. The interchange between integration and differentiation used to obtain (12 8 p is 
standard. One can justify it properly by taking the integral over a compact interval of 
the type [0, x] first and deduce an expression for this function of a similar kind of ( !29|) . 
Taking the limit x —> 00 gives the result. 



2.3 Computation of z/2 

Due to ([SD and ([7]) 1^2 (^) and z>2(/i) are related by 

POO 

v^{h) = / e^H-^''^V2{t-^''^h)dt. (31) 
Jo 

The goal of the present subsection is inverting this integral transform. 

By straightforward computations - performing a change of variables and an integra- 
tion by parts - form (I12p and (1ST]) we obtain 

u\z"') = \ r e-'^P2{s-^")s'^"ds, 
3 Jo 

and hence 

V2{h) = ^h-\f*f){h~'), (32) 
where / is the inverse Laplace transform of the function z H- u{z^^'^): 

POO 

u{z^/^) = / e-^'f{s)ds, 
Jo 

and * stands for convolution on [0, 00). 

In the following computations we denote generically by C a multiplicative constant. 
The value of this may change from line to line. These constants could be written explicitly 
in each step, but this wouldn't be much illuminating. The value of the norming constant 
in ([T3|) will be identified at the very end of the computations. 
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The function / is explicitly known, see [11]: 

/(s) = Cs-^/3e-2/(9«). 

Thus, 

/ * f{s) = C r {t{s - t))-'/'e-(2«)/(9*(^-*))f;t 
Jo 

= Cs-'/^ f {til - t))-'/'e-2/(9^*(i-*))t^t 
Jo 

POO 

Jo 

= Cs-^/3e-^/(^^)f/(l/2, 4/3; 8/(9s)). 

= Cs-^/=^e-^/(9^)f/(l/6, 2/3; 8/(9s)). (33) 

In the third line we perform the change of integration variable (t(l — t))~^ =: 4(1 + z). 
Finally, in the last step we apply Kummer's identity (153|) . 
From (152]) and (I33p we readily get 

zy2(/i) = Ce-(«'^')/V(l/6, 2/3; (8/i3)/9). 

Identification of the norming constant C in the last expression follows from (fT2|) . fj44|) 
and (IM]): 

24/33i/3r(2/3)r(l/3)~i = z>2(0) = r(2/3)z/2(0) 

= CT(2/3)f/(l/6,2/3;0) = Cr(2/3)r(l/3)/r(l/2). 

Hence 

2.4 Computation of z>i 

Recall the definition of Ui: 

/■oo 

z>i(x) = / z>(a;, h)dh. 
Jo 

Let us examine the case x > (the function Oi is even). Using (127|) we readily obtain 

1 /'°° 

(9^z>i(a;) = 2 y («(x)'9hV^(x, x) - u'{x)vix, x)) dx 

= i(u'(0)u;(x)-M(0)4v^(x,0)). (34) 
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Remark. Again, the interchange between integration and differentiation is classicaL 
Note first that for every x > 0, dh<f{x, /i) — )■ when h — i- oo: It admits a finite or infinite 
hmit (integrating (l23l) with respect to h, one can see that the terms involved admits finite 
or infinite hmits - recall dx^p{x, h) is negative) and its integral is finite. One can then 
integrate dxh'{x^h) first on [0, /i] and obtain an expression for J^^ u{x)^{x,x)dx- Taking 
/i — 7- oo and using lim/i_j,oo dh<f{x, h) = rigorously prove (134|) . 

Notice that we have dhif{x,0) = for every x (see (124|) ). The following lemma will 
serve to compute w{x) on the right hand side of (134|) . 



Lemma 1. The function (f : R_|_ x R_|_ — )■ R defined in (l^T]) is expressed as 

-1 
u{X)u'{X) 



^(^' = -TTT^7Tt\ «(A + h) I u{X + xMx)dx (35) 



+ u{X + h) vx{X + x)u{x)dx 
Jo 

POO 

+ vx{X + h) / u{X + x)u{x)dx 



where f ^ : R — ?■ R the solution of the Airy equation ( 142|) with initial conditions 

vxiX)=u{X), v',iX) = -u'iX). (36) 

Proof of Lemma [II Let f a : R — t- R be the solution of f H2|) with boundary conditions 
(IMl) and 

(j)+{X, h) := u{X + h)l{h > 0} + vx{X - h)l{h < 0}, 
0_(A,/i) :=0+(A,-/i). 

Then 0±(A, ■) are exactly the unique solutions of (Hn|) . with \^(/;.) = such that 
lim/i^±oo 0±(A, /i) = 0. Applying (HB we readily obtain (135|) . □ 



Using now (135|) . we obtain 



w{X) = ^(A, 0) = I u{X + x)u{x)dx 



-1 



A^i'(A) 7o 
-1 



Xu\X) 

u'{X)~u'{<d)u{X) 
Am'(A) 



(u"(A + x)u{x) - u{X + x)M"(x))f^X 
[u\X + x)m(x) - «(A + x)u\x))'dx 

(37) 
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In the first step we have used (1351) . in the second step the Airy equation 
Now, using the key identity (H6i) and the trace formula (H7I) . from (137|) we derive 



2 -^fc + A' 

k=l 



and hence 

|«'(0)| 



w[x) 



2 

A;=l 



(38) 



Finally, dSl), dSH]) and ([381) yield (ffll)- 



2.5 Computation of ui 

Due to ([6]) and ([7]) z^i(a;) and vi{x) are related by 

/■oo 

i)^[x) = / e-H~'^'^iJi{r'^/^x)dt. 
Jo 

The goal of the present subsection is inverting this integral transform. Since 1^1(0;) is 
convex combination of scaled exponentials, we only have to compute the inverse transform 
of the exponential density. Assume 

POO 

e-^= / e"*t"2/^0(t"2/3x)dt, a;>0. 
Jo 

Writing the integer moments of both sides, by elementary computations we obtain 



(f){x)x^dx 



r(2n/3 + l)' 

This identifies (f) as the Mittag-Leffler density /2/3 from f l5T]) and fl54p thanks to f l5U]) . 
Hence ([151). 



3 Appendix: 

3.1 Feynman-Kac formulas 

According to the present context we formulate the Feynman-Kac formulas only for the 
one- dimensional Brownian motion. For the general theory of Feynman-Kac formulas see 
e.g. [9J, [lOj, or [17J. 

Let y B{y) he a standard 1-dimensional Brownian motion which starts from level 
B{0) = h eM, and V : R — t- ]R_|_ and / : R — t- M be continuous functions. Assume that / 
is also bounded. Define 

g{x,h) := E(e-io'^(^(J'))'^V(5(a;)) | ^(0) = h), x> 0, 

POO 

g{\, h) := / e'^^'gix, h)dx, A > 0. 

Jo 
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Theorem. [The Feynman-Kac formula for 1-d Brownian motion] 

(i) {x,h) I—)- g{t, x) is the unique bounded solution of the parabolic Cauchy problem 

d,g{x, h) = Ulg{x, h) - V{h)g{x, h), g{0, h) = f{h). 

(a) For A > fixed, h i— )■ ^(A, h) is the unique bounded solution of the elliptic PDF 

^dl~g{X, h) = {V{h) + A)^(A, h) - f{h). (39) 

(Hi) Fix A > and let 0±(A, ■) : M — > M be the unique solutions of the ODF 

^dim = {v{h) + x)m, (40) 

for which which \imh^±oo 4'±{X, h) = 0. (The convergence is actually at least exponentially 
fast.) Then the unique bounded solution of (139|) is expressed as: 

~g{\, h) =2 (0+(A, O)0'_(A, 0) - 0_(A, O)0V(A, 0))"' (41) 
x|0+(A,/i) (t>-{\x)f{x)dx-ct>-{\x) £ M\x)f{x)dx^- 

For a full proof of these statements see e.g. [8]. 
3.2 Airy functions 

In the present Appendix we collect the necessary minimum information about Airy func- 
tion needed and used in this paper. For an exhaustive treatment of Airy functions see 

Consider the second order linear ODE for / : M — > M: 

f"{h) = 2hf{h), (42) 

commonly called Airy's equation. We denote by -u : M — t- M the unique solution of (142|) 
with boundary conditions 

m(0) = 1, lim u{h) = 0. (43) 

In terms of the conventionally defined Airy function of the first kind Ai(-) we have 

u{h) = 32/=^r(2/3)Ai(2^/=^/i), 
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Figure 6: The normalized Airy function u 



and we also have 



uih) ~ ^^P^h-y'exp {-i2'/y3)h'/') , as ^ oo (45) 



uih) ~ ^^^^^ \h\-'/' sin ((2^/73) \hf' + 7r/4) , as ^ -00. 



The function m : R — )■ M is convex and decreasing on [0, 00) and oscillates indefinitely on 
(—00, 0]. The function u extends as an entire function to C 3 z u{z) G C. 

We denote by k = 1,2, . . . the consecutive zeros of h u'{—h). It is known that 

0<S[<6',<---<6',<..., 
and their asymptotics, as k ^ 00, are 

5'k - ^(Svrfc)^/^ 

(Finer asymptotics are also available but not needed for our purposes in this paper.) 
The following key identity holds 

It is easily seen that the sum on the right hand side is absolutely convergent and the two 
sides have simple poles at —5^, k = 1,2, . . . , with the same residues. Thus, the two sides 
can differ in an entire function only. For a full proof of this formula see [6j. 

From (H6|) an infinite family of trace formulas for YlT=i(^'k)''"^ — 2) are expressed 
in terms of u'{0). We will use only the first three of these: 

00 „ 00 00 „ 

&ir'-% (47) 

k=l ^ ' k=l k=l ^ ' 
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These are easily checked by differentiating both sides of (1471) at z = and using 
We will denote 

■■= ^(S'k)-'- (48) 

Note that Yj'k=iPk = 1- 

3.3 Mittag-Leffler distributions 

Following the terminology of [5j and [3j, we will call Mittag-Leffler distribution of index 
a G [0, 1] the probability distribution : M_|_ — )■ [0, 1] whose moment generating function 
(Laplace transform) is 

e-y^dF^{x) = Y: fti^, =■■ Eai-y). (49) 



k=0 



The function y i-> Ea{y) defined by the power series in (H^ is the Mittag-Leffler function 
of index a G [0, 1]. It is proved in [15j - where reference is also made to unpublished 
alternative proof of W. Feller - that, for a G [0,1], [0, cxd) 3 y v-^ Ea{—y) is indeed 
completely monotone, and thus the Laplace transform of a probability density function 
on R+. 

The moments of the Mittag-Leffler distributions are 

x"^dF^{x) = — -. (50) 

r(am + l) 



a = 0, |, 1 are special cases: 

Fo{x) = 1 - e-^ Fi/2(x) = \/| /" e-y"/^dy, Fi{x) = l{x > 1}. 



For a G [0, 1) the distribution function F^ is smooth and the following power series 
expansion holds for the density function /^(x) := F^{x), see 



Ux) = -y ^'^((k + l)an)m + l)a)^_^^, ^^^^ 

k=0 



This power series defines actually an entire function on C. See the picture below for a 
representation of those functions with a = 0, 1/2 and 2/3. 

We are interested in the case a = 2/3. This particular case is also expressed in terms of 
a confluent hypergeometric function, see (15^ below. 
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Figure 7: Mittag-Leffler density functions for a 
a = 2/?> (thick) 



(dotted), a = 1/2 (plain) and 



3.4 Confluent hypergeometric functions 

See [7j, [IJ for details on confluent hypergeometric functions. We collect here only a few 
facts needed for our purposes. 

Let a > and 6 G M be flxed parameters. We deflne the confluent hypergeometric 
functions of the second kind, or Tricomi functions 



U (a, 6; 



by the integrals 



U (a, h] z) 



(52) 



r(a) JO 

(For more general deflnitions see [7], jl].) In our formulas the particular cases t/(l/6, 2/3; ■), 
f/(l/2,4/3; ■), and f/(l/6,4/3; ■) will occur. 
Kummer's identity 

U{a,b;z) = z^-''U{l + a-b,2-b,z) (53) 

holds if 6 > a + 1. 

The Mittag-Leffler density function /2/3, which appears in (fT5|) is also expressed in 
terms of f/(l/6,4/3; ■) as follows: 

f2/s{z) = ^ze-^'^->^'''U{l/6,4/S; (4z^)/27). (54) 



)/2'?7(l/6,4/3; (4z3)/27). 

Computing the moments of the right hand side of flM|) we obtain the expressions (jH)]) . 
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